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ABSTRACT 

\ We use high-resolution cosmological hydro dynamical AMR simulations to predict the 

characteristics of La emission from the cold gas streams that fed galaxies in massive 
haloes at high rcdshift. The La luminosity in our simulations is powered by the 
, release of gravitational energy as gas flows from the intcrgalactic medium into the halo 

Q_i' potential wells. The UV background contributes only < 20% to the gas heating. The 

JL I La emissivity is due primarily to electron-impact excitation cooling radiation in gas ~ 

5_i . 2 x 10 4 K. We calculate the La emissivities assuming collisional ionisation equilibrium 

(CIE) at all gas temperatures. The simulated streams arc self-shielded against the UV 
5j ! background, so photoionisation and recombination contribute negligibly to the La line 

formation. We produce theoretical maps of the La surface brightnesses, assuming that 
~ 85% of the La photons are directly observable. We do not consider transfer of the La 
radiation, nor do we include the possible effects of internal sources of photoionisation 
such as star-forming regions. Dust absorption is expected to obscure a small fraction 
I of the luminosity in the streams. We find that typical haloes of mass M v ~ 1O 12_13 M0 

, at z ~ 3 emit as La blobs (LABs) with luminosities 10 43-44 ergs -1 . Most of the La 

■ comes from the extended (50—100 kpc) narrow, partly clumpy, inflowing, cold streams 

of (1 — 5) x 10 4 K that feed the growing galaxies. The predicted LAB morphology is 
therefore irregular, with dense clumps and elongated extensions. The integrated area 
' contained within surface-brightness isophotcs of 2 x 10~ 18 ergs" 1 cm -2 arcscc -2 is 

f*^ \ ~ 2 — lOOarcsec 2 , consistent with observations. The linewidth is expected to range 

from 10 2 to more than 10 3 kms~ 1 with a large variance. The typical La surface 
brightness profile is oc r -12 where r is the distance from the halo centre. Our simulated 
LABs are similar in luminosity, morphology and extent to the observed LABs, with 
jJJ ■ distinct kinematic features. The predicted La luminosity function is consistent with 

observations, and the predicted areas and linewidths roughly recover the observed 
scaling relations. This mechanism for producing LABs appears inevitable in many 
high-z galaxies, though it may work in parallel with other mechanisms. Some of the 
LABs may thus be regarded as direct detections of the cold streams that drove galaxy 
evolution at high z. 

Key words: cosmology — galaxies: evolution — galaxies: formation — galaxies: high 
rcdshift — intergalactic medium — galaxies: ISM 
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Hundreds of Lyman-alpha blobs (LAB) have been detected 
so far in the redshi f t range z = 2 — 6.5, mostly near z ~ 3 
JSteidel et al.ll200d : iMatsuda et alj|2004 ISaito et alj 120061 : 
lOuchi et al.l 120091 : lYang et al.l l20od h Their luminosities 



© 2009 RAS 



2 Goerdt, Dekel, Sternberg, Ceverino, Teyssier & Primack 



range from below 10 43 to above 10 44 ergs _1 , and they 
extend on the sky to 30 — 50kpc and more. The two main 
open questions are (a) the origin of the extended cold and 
relatively dense gas capable of emitting La , and (b) the 
continuous energy source for exciting the gas to emit La . 

The emitting hydrogen gas should be at a temperature 
T Si 10 4 K, relatively dense and span a much larger area than 
covered by the stellar component of galaxies. It may arise 
from outflows or inflows. The energy sources discussed in the 
literature include photoionisation by obscu red AGNs, early 
starbursts or extended X-ray emission dHaiman fc Reesl 
120011 : 1 Jimenez fe Haimanl [2009 ; IScharf et all 120031 ). as well 
as compression of ambient gas by superwinds to dense 
La emitting shells |Mori et al.l |2004 ) , a nd star for mation 
triggered by relativistic jets from AGNs (I Reesl 1 19891 ) . 

Many of the bright LABs are found in the vicinit y 
of massive, star- forming galaxies (|Matsuda et al.l 120061 ). 
Multi-wavelength observations reveal that a fraction of 
the LABs are associated with sub-millimetre and infrared 
sources that indicate very high star- formation rates (SFR 
in the ranee 10 2 - lO 3 Af yr _1 (|Chapman et al.l 



(SFR) 
1 l200ll ; 



a.ng e 

iGeach et a.1.1 120051. 2007|) or with ob s cured active galactic 
nucle i (AGN) (|Basu-Zvch fc ScharJ 12004 IScarlata etail 



120091 ). While stellar feedback and AGNs could in principle 
provide the energy source for the La luminosity, many 
LABs are not associated with sources of this sort that are 
powerful enough to explain the observed La luminosities 
(|Nilson et all I2OO6I ; ISmith fc Jarvisl Il999l ). This indicates 
that star formation and AGNs are not the sole drivers of 
LABs, and may not even be the dominant source for LABs. 

Indeed, high-redshift 
galaxies exhibit a generic mechanism that simultaneously 
provides both the cold gas and the energy source. It is 
a direct result of the phenomenon robustly established by 
simulations and theoretical analysis, where high-z massive 
galaxies are continuously fed by narrow, cold, intense, 
partly clumpy, gaseous streams that penetrate through the 
shock-heated halo gas into the inner galaxy, grow a dense, 
unstable, turbu lent disc with a bulge and trigger rapid 
star formation (iBirnboim fc Dekel 20031; iKeres et al.l 120051; 



Keres et al.l 



20091; 



Dekel fc Birnboiml 120061; lOcvirk. Pichon fc Tevssierl 



Dekel et al 



2008 



2009 



Johansson. Naab fc Qstrikerl 200E ; Dekel. Sari fc Ceverinol 



20091 ; ICevermo. Dekel fc Bournaudlf2009jV_ Massive c l ump; 
star-forming disk s observed at z ~ 2 (IGenzel et alj I200S., 
iGenel et afTl2008l ; iForster Schreiber et al.l 120091 1 may have 
been formed via the smooth and steady accretion provided 
by cold flows, as opposed to merger events. The streaming of 
the gas into the dark-matter halo potential well is associated 
with transfer of gravitational energy to excitations of 
the hydrogen atoms followed by c o oling emission of 
La dHaiman. Spaans fc Quat acrt 2000£_ Fardal et al.l l200ll 



Dckcl & Birnboiml 20061. boOSl ; iKhochfar fc Ostrikeil [200 



Furlanetto et al J 12003. 2005). 



There were two earlier attempts to compute the La 
cooling radiation from hydrodynami c SPH cosm ol ogica l 
simulations. Based on their analysis, iFardal et all (|200ll ) 
pointed out the potential association of this La cooling 
emission with the first observed LABs. These early 
simulations did not allow a proper resolution of the detailed 
structure of the cold streams as they penetrate through 
the hot medium. Their shortcomings included intrinsic 



limitations of the SPH technique, a limited force resolution 
of 7/i _1 kpc (comoving), not allowing for radiative cooling 
below 10 4 K and neglecting the pho toionisation by the UV 
background. iFurlanetto et al] (|2005l ) used SPH simulations 
with a somewhat higher resolution of ~lkpc to make 
more detailed comparisons of the simulated and observed 
luminosity functions and size distributions of LABs powered 
by cold accretion. In their pessimistic model, assuming no 
cmissivity from gas that is self-shielded from the ionising 
background radiation, they end up with low La luminosity. 
They comment that cooling IGM gas may explain the 
observations only if one adopts an optimistic scenario, where 
the self-shielded gas is emitting La at CIE. The latter 
scenario will also be adopted in this paper. Other sources of 
ionisation are ignored and radiative transfer is not applied. 

iDiikstra fc Loebl (|2009l . hereafter DL09) have recently 
worked out an analytic toy model for La cooling radiation 
from cold streams, based on the general properties of the 
cold streams as reported from cosmological simulations. 
They conclude that the streams could in principle provide 
spatially extended La sources with luminosities, line widths 
and abundances that are similar to those of observed 
LABs. They point out that the filamentary structure 
of cold flows may explain the wide range of observed 
LAB morphologies. They also highlight the fact that the 
most luminous cold flows are associated with massive 
haloes, which preferentially reside in dense large-scale 
surroundings, in agreement with the observed presence of 
bright LABs in dense environments. This model presents 
a successful feasibility test for the role of cold streams 
in powering the LAB emission, and it provides physical 
intuition into the way by which this process could be 
manifested. However, a comparison to our simulations 
indicates that this simplified model does not capture the 
detailed hydrodynamic properties of the cold streams. In 
particular, it significantly over-predicts the cold-gas density, 
under-predicts the gas temperature, and does not account 
for the partly clumpy nature of the streams and their 
characteristic radial distribution in the halo. As a result, 
the DL09 model underestimates the total La luminosity by a 
factor of a few, and it therefore has to appeal to the excessive 
clustering of the LABs in order to boost up the predicted 
luminosity function for a match with the observations. 

In the current paper, we calculate the La emission 
from the cold gas in high-redshift galactic haloes 
using state-of-the-art hydrodynamical AMR cosmological 
simulations. With a maximum resolution better than 70 pc 
in our simulations we can quite accurately map the extended 
La sources. This allows us to study their individual shapes, 
morphologies and kinematics. We measure quantities such as 
the distribution of surface brightness within each halo, the 
area covered with surface brightness above a given isophotal 
threshold, the predicted observed linewidth, the typical total 
La luminosity per halo of a given mass, and the overall La 
luminosity function of LABs. These predicted properties are 
compared with the observed LABs. 

This paper is organised as follows. In fJ3] we work 
out a simple feasibility test where we use a simple toy 
model to estimate the expected gravitational heating power 
as a proxy for the La luminosity. In iJS] we introduce 
the two sets of cosmological simulations used. In ||4] 
we explain our methodology of computing La emissivity 
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and luminosity as a function of halo mass. In ij5] we 
identify the gas that contributes to the La emission. In 
Sj6] we apply our methodology to the simulations and 
provide predicted images of LABs. In iJ7]we determine the 
scaling of La luminosity as a function of halo mass and 
redshift. In ij8] we compare our predicted La luminosity 
function with observational results. In ij9] we measure 
the predicted surface-density profile, isophotal area and 
linewidth and compare to observations. In H 101 we show that 
gravitational heating is the main source of energy driving 
the La luminosity in our simulations, while the role of 
photoionisation is minor. In Ulll we discuss our analysis and 
results and draw our conclusions. 



2 FEASIBILITY OF GRAVITATIONAL 
HEATING 

The gravitational energy gain due to the streaming of the gas 
from the virial radius toward the centre of the halo potential 
well is a natural source of energy for the La emission. The 
gravitational energy released per unit infalling mass is of 
order V 2 = GM V /R V , where the quantities are the virial 
velocity, mass and radius respectively. The accretion rate 
M can be estimated from the observed star-formation rate. 
For V ~ 300kms _1 and M ~ 150M Q yr _1 we obtain a 
power M V 2 ~ 10 43 ergs _1 , comparable to the luminosity 
of a typical LAB. 

In some more detail, an analysis of the MN cosmological 
simulation (described below in ^3} reveals that both the cold 
gas accretion rate M c and it s inward velocity are roughly 
constant along the streams (|Dekel et al.l |20(M ). If so, the 
gravitational power deposited at radius r per unit radial 
length in the cold gas is 



-Egrav — f c M c 



dr 



(1) 



where cj>{r) is the gravitational potential at r. The factor / c 
is the fraction of the energy that goes to heating the cold 
streams themselves rather than the hot medium. The total 
power deposited between the virial radius R v and radius r, 
across a potential difference A<f>(R v ,r), is 



f c M c \A^(R v ,r)\. 



C2) 



Given the total mass density profile p(r) and the 
associated mass profile M(r), the potential gain is 



-1 + 



V 2 (r) 
V 2 



Mr') 



r'dr', 



(3) 



where V 2 (r) = GM(r)/r, p v is the mean density within the 
virial radius, and (j> is typically a number of order unity. 

For an NFW potential well with a concentration 
parameter C (|Navarro. Frenk fc Whit3ll997h . one has 



(r) 



C 



MC) 



ln(l + x) ln(l + C) 



x 



c 



(4) 



where x = Cr/R v and Ai(x) = ln(x + 1) — 
x/(x + 1). Accordin g to cosmological TV-body simulations 
|Bullock et al.ll200ll h the average concentration parameter 
as a function of halo mass and redshift is C ~ 3Mf 2 0J3 (l + 
z)7\ where M i2 = M V /W 12 M Q and (1 + z) 4 = (1 + z)/4. 




r/R„ 



Figure 1. Cumulative luminosity profile, showing the fraction of 
the luminosity from within a sphere of radius r compared to the 
virial sphere. The toy-model prediction (solid red) is proportional 
to the fraction of gravitational energy that is deposited inside 
radius r based on eq. (JSJ. The predicted profile in the range (0.2 — 
0.5)i? v can be approximated by the power law L(< r) oc r 0,7 . 
About half the luminosity is expected to originate from the cold 
streams in the halo outside the inner 0.27R V sphere. Shown in 
comparison is the luminosity profile as derived from the three 
simulated CDB galaxies stacked together (see section [T0)| . The 
similarity between the toy model predictions and the simulation 
results is remarkable. 



For M v = lO 12 A/ , this is C ~ 3 at z = 3. Then <f>(r) ~ 
4.7 [ln(l + x)/x- 0.46] so $(x -> 0) ~ 2.5. 

A practical approximation for the average accretion rate 
into haloes of M v ~ 10 12 ~ 13 Af o in the standard ACDM 
cosmologjQ is derived from the EPS approximation and from 

the Millennium cosmologica l simulation 

( Neistein. van den Bosch fc Dekel 120061; iNeistein fc Dekell 



l200Sl : Birnboim, Dckcl fc Neist ein 2007j, Appendix A). The 
same MN hydrodynamical simulation that is used in the 
current paper (33J) reveals that 95% of the gas accretion is 
through cold streams that penetrate efficiently all the way 
into the vicinity of the central galaxy. For a baryonic fraction 
of 0.165 in the incoming streams, the average cold accretion 
rate is approximated by 



137 Mqyt- 1 Ml 2 15 (1 + z)f 25 . 



(5) 



This accretion rate is consistent with the typical SFR of 
~ lOOA/Qyr -1 observed in massive galaxi es of similar 
comoving number densities at redshifts 2 — 3 dGenzel et al.l 
20061; iFoister Schreiber et al.l l200d: lElmegreen et al.ll2007l ; 



Genzel et al.ll2008l ; IStark et al.ll20oj h 

The virial velocity is given b>0 



1 Adopting the parameters motivated by WMAP5: Q n 
fl A = 0.72, h = 0.7, cr 8 = 0.8. 

2 With R v ~ 77kpcMj^ 3 (1 + z)^ 1 , for completeness. 



0.28, 
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V v ~ 236 kins" 1 M 1 / 3 (1 + z) 
so we finally obtain in eq. ([2| 



1/2 



E a 



1.2 x 10 43 ergs -1 f c M^ 82 (1 + z) 



3.25 
4 



(6) 



(7) 



Assuming that a substantial fraction f a of this energy is 
emitted as observable La radiation, we conclude that at 
z ~ 3, with f a fc ~ 1, a LAB of luminosity 1/ ~ 10 43 ergs -1 
is feasible from haloes of M v ~ 10 12 Mo. Luminosities as 



10 1 

high as ~ 10 ergs" require haloes of ~ 3 x 1O 12 M . 
If fafc is only ~ 0.1, then the required haloes for the 
same luminosities are about three times more massive. 
We note that the mean comoving number density for 
haloes more massive than (1,3,10) x 1O 12 M0 at z — 3 
is (4.7,0.68,0.057) x 10~ 4 Mpc~ 3 respectively. Given that 
LABs of a luminosity higher than 10 43 ergs -1 appear with 
a comoving number density of ~ 5 x 10~ 5 Mpc~ 3 (Fig. I14[) . 
the simple gravitational heating model indicates that they 
can emerge from haloes of ~ 3 x 1O 12 M0, with f a fc ~ 0.1. 
Thus, the comparison of the estimates from our toy model 
with the total luminosities of the observed LABs indicates 
that gravitational heating is a feasible source for the La 
emission. This kind of energy has to be released from these 
galaxies. 

One can combine equations @ and (Q) to evaluate the 
gravitational energy deposited at different radii, 



1 / c A/ 1 2 82 (l + 2) 3 - 25 



£grav(<r) ~ 1.2 x 10 43 ergs 

1.861n(l + Cr/-R v ) 

86 c7JW v 



(8) 



This gives a very crude estimate for the 3D luminosity profile 
L(< r) — / a /ci?grav(< r). This profile is shown in Fig. [TJ 
normalised to the total luminosity inside the virial radius. 
One can read from this plot the fraction of the luminosity 
that is expected in the different zones of the halo. For 
example, about half the luminosity is expected from outside 
r = 0.27R V , implying that in haloes of M v ~ 10 12 M Q the 
LABs are expected to extend over more than 50kpc. The 
luminosity profile predicted by gravitational heating is to 
be compared with our results from the simulations (S UOI) . 
and with observed LABs (3HJ). 



3 SIMULATIONS 

We use here simulated galaxies from two different suites of 
simulations, both employing Eulerian Adaptive Refinement 
Tree (AMR) hydrodynamics in a cosmological setting. 
One suite consists of three simulations zooming in with a 
maximum resolution of 35 — 70 pc on individual galaxies 
that reside in dark-matter haloes of masses ~ 5 x 1O 11 M0 
at z = 2.3 (jCeverino. Dekel fc Bournaudl 120091 . hereafter 
CDB ) . The other suite is from the Horizon MareNostrum 
simulation containing hundreds of massive galaxies in a 
cosmological box of sid e 50fe -1 Mpc with a maximum 
resolution of ~ 1 kpc (|Ocvirk. Pichon fc Tevssierl 120081 . 
hereafter MN). 

Figure [2] shows sample gas density maps of galaxies 
from the two suites of simulations. They demonstrate 
the dominance of typically three, narrow cold streams, 
which come from well outside the virial radius along the 
dark-matter filaments of the cosmic web, and penetrate 



into the discs at the halo centres. The streams are partly 
clumpy and partly smooth, even in the simulation of higher 
resolution. The typical densities in the streams are in the 
range n = 0.01 — 0.1 cm -3 , and they reach n = 0.1 — 1 cm -3 
at the clump centres and in the central disk. 



3.1 High-resolution ART CDB simulations 

The CDB simulations have been run with the AMR tree 
code Art l|Kravtsov. Klypin fc Khokhlovl 1 19971 ; iKravtsovl 
l2003h with a spatial resolution better than 70 pc in 
physical units. The code incorporates the relevant physical 
processes for galaxy formation, including gas cooling and 
photoionisation heatin g, star formation, metal e nrichment 
and stellar feedback (|Ceverino fc Klypinl l2009t ). Cooling 
rates were computed for the given gas density, temperature, 
metallicity, and U V background based on Cloudy 
(|Ferland et al.l Il998l ). assuming cooling at the centre of 
a cloud of thickness 1 k pc (|Ceverino- Rodriguez! 120081 ; 
iRobertson fc Kravtsovl2008l ). Metals are included, assuming 
an abundance of 0.02 relative to the total mass for a solar 
composition. The code implements a "constant" feedback 
model, in which the combined energy from stellar winds and 
supernova explosions is released as a constant heating rate 
over 40 Myr, the typical age of the lightest star that explodes 
as a type-II supernova. Photo-heating is also taken into 
account self-consistently with ra diative cooling. A u niform 
UV background based on the lHaardt fc Madaul (|l996l ) 
model is assumed, local sources are ignored. In order to 
mimic the self-shielding of dense, galactic neutral hydrogen 
from the cosmological UV background, the simulation 
assumes a substantially suppressed UV background (5.9 x 
10 26 erg s -1 cm~ 2 Hz -1 , the value of the pre-reionisation 
UV background at z = 8) for the gas at total gas densities 
above n = 0.1 cm -3 . 

The unique feature of this code for the purpose 
of simulating the detailed structure of the streams and 
the gravitational instability in the disk is to allow the 
gas to cool down well below 10 4 K, thus reaching high 
densities in pressure equilibrium with the hotter and 
more dilute medium. A non-thermal pressure floor has 
been implemented to ensure that the the Jeans length 
is resolved by at least seven resolution elements and 
thus p revent artificial fragmentation on the smallest grid 
scale dTruelove et alJ Il997l: IRobertson fc Kravtsovl [2008; 
ICeverino. Dekel fc Bournaudl 120091 ). It is effective in the 



dense (n > 10 cm" ~ 3 ) and cold (T < 10 4 ^) regions inside 
galactic disks, while most of the La emitting gas, which is at 
temperatures (1 — 5) x 10 4 K, is not affected by this pressure 
floor. 

The equation of state remains unchanged at all 
densities. Stars form in cells where the gas temperature 
is below 10 4 K and its density is above the threshold n = 
lcm -3 according to a stochas ti c mo del that is roughly 
consistent with the iKennicuttl (| 19981 ) law. The ISM is 
enriched by metals from supernovae type II and type la. 
Metals are assumed to be released from each stellar particle 
by SNII at a constant ra t e for 40 Myr since its birth, 
assuming a iMiller fc Scald (Il979h IMF and matching the 
results of lWooslev fc Weaver! ( 19951 ). The metal ejection by 
SNIa assumes an exponentially declining SNIa rate from a 
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Figure 2. Gas density in simulated galaxies from CDB (left) and MN (right). The colour refers to the maximum density along the 
line of sight. The contours mark n = 0.1, 0.01 and 0.001 cm~ 3 , respectively. The circle shows the virial radius. Left: one of the three 
CDB galaxies (resolution 70pc) at z = 2.3, with M v = 3.5 X 10 11 Mq. Right: one of the MN galaxies (resolution 1 kpc) at z = 2.5, with 
jtf v = 10 12 M Q . In both cases, the inflow is dominated by three cold narrow streams that arc partly clumpy. The density in the streams 
is n = 0.003 — 0.1 cm -3 , with the clump cores reaching n ~ 1 cm -3 . 



self-consistently (|Ceverino-Rodriguezll2008l ) . 

The dark matter particle mass is 5.5 x 10 M®, the 
minimum star particle mass is 10 Mq, the smallest cell size 
is 70 pc (physical units), and the force softening length is 
105 pc. 

The initial conditions for the CDB simulations 
were created using a low-resolution cosmological iV-body 
simulation in a comoving box of side 20/i~ 1 Mpc, for which 
the cosmological parameters were motiv ated by the WMAP5 
following values (|Komatsu et al.ll2009l ): Sl m = 0.27, SI a = 
0.73, Q b = 0.045, h = 0.7 and erg = 0.82. At z = 1, three 
haloes of M v ~ 10 12 M Q each have been selected, avoiding 
haloes that were subject to a major merger near that time. 
The three halo masses at z = 2.3 are 3.5,4,6 x 10 Mq, 
and they end up as (3 — 4) x 10 12 M© haloes today. For each 
halo, a concentric sphere of radius twice the virial radius was 
identified for re-simulation with high resolution. Gas was 
added to the box following the dark matter distribution with 
a fraction /b = 0.15. The whole box was then re-simulated, 
with refined resolution only in the selected volume about the 
respective galaxy. 



maximum at 1 G yr. The code treats the ad vection of metals photo- ionisation (|Katz. Hernquist fc W einberg 1992). Metal 

cooling is also included using tabulated Cloudy rates, 
and is assumed proporti o nal t o the metallicity, relative 
to the iGrevesse fc Sauva 3 (|l99Sh solar abundances. Unlike 
in the CDB simulation, no cooling below T < 10 K is 
computed, and no self-shielding of the UV flux is assumed. 
Because we assume that the Lev emissions are produced 
in predominantly shielded gas that is in CIE, our MN 
computations should be regarded as less accurate than the 
CDB simulation. 



3.2 RAMSES Mare Nostrum simulations 

The MN sim ulation uses the AMR code Ramses 
jTevssierl l2002t ). The spatial resolution is ~ lkpc in 
physical units. UV he ating is included assuming the 
lHaardt fc Madaul (1 19961 ) background model, as in the 
CDB simulation. The code incorporates a simple model 
of supernovae feedback and metal enrichment u si ng th e 
implementation described in iDubois fc Tevssierl (|2008l ). 
The cooling rates are calculated assuming ionisation 

and 



For high-density regions, the MN code considers 
a polytropic equation of state with 70 = 5/3 to 
model the complex, multi-phase and turbulent structure 
of the inter-stellar mediu m (ISM) (|Yepes et al ] Il997l : 
Springel fc Hernquis tl 1 19991) in a simplified form (see 



Schave fc Dalla Vecchiall999l ; lDubois fc Tevssierl2008l ). The 



ISM is defined as gas with hydrogen density greater than 
71h = 0.1 cm' 3 , one order of magnitude lower than in the 
CDB simulation. Star formation has been included, for ISM 
gas only, by spawnin g star particles at a rate consistent with 
the lKennicuttl (| 19981 ) law derived from local observations of 
star forming galaxies. 

The MN simulation implemented a pressure floor in 
order to prevent artificial fragmentation, by keeping the 
Jeans lengthscale, Aj oc Tn~ 2//3 , larger than 4 grid-cell 
sides everywhere. In any case where n > 0.1cm -3 , 
a density dependent temperature floor was imposed. It 
mimics the average thermal a nd turbulent pressure of the 
m ultiphase ISM, in the spir i t of Springel fc Hernouistl [1 19991 ) 



m ultiphase Iblvl, m the spir it or ap: 
or lDalla Vecchia fc Schavel (|200ct ). 



are calculated assuming 
equilibrium for H and He, including both collisional- 



In our case, we allow the 
gas to heat up above this temperature floor and cool back. 
The temperature floor follows a polytropic equation of state 
with T floor = T (n/no)' lo ~ 1 , where T = 10 4 K and n = 
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0.1 atoms cm 3 . The resulting pressure floor is given by 

-Pfloor = ^HfcBlfloor- 

We crudely correct for this artificial temperature 
boost in post-processing by subtracting Tfl oor from the 
temperature read from the grid cells where n > 0.1cm -3 . 
If the corrected temperature is below 10 4 K, we set it to 
10 4 K. In practise, almost all cells where n > 0.1cm -3 are 
set to T = 10 4 K. As will become clear in fjll this means 
neglecting any La emission from these high-density cells. 
We will evaluate the possible error made by this procedure 
in the MN galaxies using the more accurate high-resolution 
CDB galaxies, where no temperature floor has been applied. 

For each stellar population, 10% of the mass is assumed 
to turn into supernovae type II after 10 Myr, where the 
energy and metals are released in an impulse. For each 
supernova, 10% of the ejected mass is assumed to be pure 
metals, with the remaining 90% keeping the metallicity of 
the star at birth. SNIa feedback has not been considered. 

The dark matter particle mass is 1.16 x 10 7 Mq, the star 
particle mass s 2.05 x 1O 6 M0, the smallest cell size is 1.09 
kpc physical, and the force softening length is 1.65 kpc. 

The initial conditions of the MN simulation were 
constructed assuming a ACDM universe with Hm = 0.3, 
n A = 0.7, fib = 0.045, h = 0.7 and erg = 0.9 in a 
periodic box of 50 /i -1 Mpc. The adaptive-resolution rules in 
this simulation were the same everywhere, with no zoom-in 
resimulation of individual galaxies. 



4 COMPUTING THE La LUMINOSITY 

Given the gas temperature T and mass density p in 
every cubic grid cell of the hydrodynamic simulation, we 
compute the local La emissivity produced by electron 
impact excitation of neutral hydrogen, 



e = n c n m qi s ^2 P (T) huL a erg cm s 



(9) 



In this expression, n e and nm are the local electron and 
atomic hydrogen particle densities (in cm -3 ) and his^a — 
10.2 eV = 1.63 x 10 -11 erg is the La photon energy. The 
temperature dependent qu antity gis-»2 P is the collisional 
excit ation rate coefficient (|Callawav. Unnikrishnan fc Ozal 
[1983), 



?ls->2p 



2.41 x 10" 



y0.5 



Go 4 ) 



x exp 



/ hULa \ 



(10) 



where T is in degrees K. Radiative recombinations of 
electrons with protons contribute negligibly to the La 
emissivity, as can be seen from the green, long-dashed curve 
in Figure [3] 

We assume a primordial helium mass fraction Y = 
0.24, corresponding to a helium particle abundance of 1/12 
relative to hydrogen. This gives 



p = (4/3) m H n H , 



(11) 



where ran is the hydrogen mass and nn is the total density 
of hydrogen nuclei (neutral plus ionised). The electron and 
atomic hydrogen densities in eq. ([9]) are given by 



WHi/fiH = xm 



(12) 



and 

n e /n H = zhii + (1/12) (x H eii + 2a- Hc 



(13) 



where we adopt the temperature dependent ionisation 
fractions, xm, imi, xnaii and XHeim as computed for 
collisional ionisation equilibrium (CIE) assuming case-B 
hydroge n recombination (O. Gn at, private communication; 
see also iGnat fc Sternben3l2007l ). A large fraction (~ 95%) 
of the La emission is produced in gas where the helium is 
fully neutral. In this limit, 



Xm = a B /(a B + qi) , 



(14) 



where ob and qi are th e recombination and co llisional 
ionisation rate coefficients ()Gnat fc Sternbe"rgll2007f l. 

We note for completeness that approximate fits can be 
provided by 



Qb (T) = 4.9 x 10"T -1& 1 



1 1 K \ -2.24 
6 m-1.5 / i , 110 \ 3-1 



yo.41 



cm' s -1 (15) 



and 

= 21.11 cm 3 s -i K 3 / 2 T- 3 / 2 exp [~) 
(2Thi/T) -1 ' 089 



[l + (5.65 T m /T) 



0.8741 1-101 ' 



(16) 



1.58 x 10 J K being the ionisation threshold 



with 1 hi = 1.5S x 1U K being the 10m 
according to lHui fc Gnedinl (119971 ). We do not use these fits 
in our calculation here. 

We assume that the cold streams are practically 
self-shielded against the photoionising UV background, for 
the following reason. The typical HI column densities along 



the shortest dimensions of the streams are 



10- 



and the total hydrogen particle densities are in the range 
0.01 - 0.1cm -3 , or typically ~ 0.03 cm -3 (see Fig. E}. At 
redshift z — 3, the mean Lyman continuu m photon intensity 
is 47r J* w 2.2 x 10 5 photons s -1 cm -2 (|Haardt fc Madaul 
Il996h and the unattenuated hydrogen photoionisation rate 
is T = 5.6 x 10 -13 s -1 . Thus, for gas densities n < 2r/a s ~ 
4.3 cm -3 , the gas will be more than 50% ionised by the 
unattenuated field. However, absorption of the radiation 
field in the stream gas will produce a photoionised column 



A, 



photo 
II 



2ttJ* 



4.2 x 10 1 ' _ 2 
cm 



(17) 

riH Ob w-h 
For the typical stream volume densities this is small 
compared to the neutral columns of ~ 10 20 cm -2 that 
we find, so the streams can be practically assumed to 
be self-shielded against the background radiation. Parts of 
the streams may not be self-shielded against local sources 
of radiation within the streams, e. g. from star-forming 
satellites. In those case we might overpredict the HI fraction 
and thus the La luminosity. 

The CIE neutral hydrogen fraction as a function of 
temperature and the temperature dependence of the La 
emissivity for a given total density are shown in Fig. [3] One 
can see that the gas is neutral at T — 10 4 K and below, and 
it becomes highly ionised at T > 2 x 10 4 K. The maximum 
emissivity at a given density is obtained at T ~ 1.8 X 10 4 K, 
and the emissivity is high enough to substantially contribute 
to the overall luminosity in the range T = (1 — 5) X 10 4 K. 
Hydrogen in this temperature range is thus the expected 
source of La emission. 
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Figure 3. Temperature dependence of the ionisation state 
and the La emissivity. The neutral-hydrogen fraction xjji 
(short-dashed, blue curve) is derived from eq. d 1 41) assuming CIE. 
The temperature dependence of the emissivity e (solid, red curve) 
is computed from cqs. J9}, 11121 and J 1 3 1 ) for n = 1 cm -3 . We see 
that the gas is neutral for T ^ 10 4 K and it becomes highly ionised 
at T > 2 X 10 4 K. The maximum emissivity at a given density is 
obtained at T ~ 1.8 X 10 4 K, and the emissivity is high enough 
to substantially contribute to an overall luminosity in the range 
T = (1 — 5) X 10 4 K. The contribution of recombination to the 
emissivity, e Q B (long-dashed, green curve), is negligible in the 
relevant temperature range. 
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Figure 4. Average metallicity profile in the simulated CDB 
galaxies, once weighted by the density of the Let emitting cold 
gas in the temperature range (1 — 5) X 10 4 K and once weighted 
by the luminosity. The vertical line indicates 0.2 R v which we 
refer to as the outer border of the galaxy. 



The total La luminosity from a given volume is assumed 
to be a constant fraction of the sum over all the cells that 
are contained in that volume, 

Ll« = fa^uVi , (18) 

i 

where Vi is the cell volume. The complex radiation transfer 
process along the cosmological line of sight is summed up in 
the transmission factor f a , the fraction of the La photons 
that make it to the observer. In i|H we will determine 
the actual value of f a by matching the predictions from 
the simulations with a preliminary determination of the 
observed La luminosity function. We expect f a to be in 
the range 0.5-1. This is because the main source of opacity 
is likely to be intervening and lower redshift intergalactic 
HI, which is expected to absorb part of the blue side of 
the line profile. An est imate for a typical line of sig ht to 
z ~ 3 is f a ~ 0.85 (e.g. iFaucher-Giguere et al.ll2008b[ ). The 
transmission factor along the line of sight to a typical LAB 
could be slightly smaller because the LABs tend to reside in 
overdense environments. A smooth component of HI in the 
emitting halo may absorb some of the red wing as well. On 
the other hand, as estimated next, dust opacity is expected 
to reduce f a by a small factor only except in the galaxy 
itself and in its immediate vicinity. 

Given the gas metallicity in every cell of the simulation, 
we can estimate the dust absorption as follows. For a 
medium in which the dust abundance scales linearly with 
the metallicity Z (with Z f» 1 corresponding to Galactic 
dust), the continuum opacity due to dust absorption at 
the La wavelength is Td us t — 0.1 Z N20 where A^o is the 
neutral- hydro gen column density in units of 10 20 cm -2 , 
(|Drainell200"3l). The line optical depth is r La ~ 10 7 N 20 T~ 1/2 
iNeufeldl (|l990h so that T du st ^ 10~ 8 Z T 4 1/2 r La . For La 
scattering through a plane-parallel medium the total line 
optical depth traversed is ~ (9aro/y / 7r) 1 ^ 3 ro, where to is 
the La optical depth from the mid-plane to the surface, 
and a = 4.7 x 10~ 4 T 4 is the damping constant (|Adams! 
Il975h . Thus, as the photon scatters through the medium the 
effective dust opacity is 

r dust ~ 2.9 Z A^ 3 T~ 1/3 . (19) 

For Z — 0.1, assuming T4 ~ 1, we obtain Td us t = 1 
for N20 — 3. We find in the CDB galaxies that the La 
luminosity-weighted fraction of volume elements that have 
a column density N20 ^ 3 between them and the observer, 
averaged over different directions, is less than a third. This 
provides an estimate of no more than a third reduction 
outside r ~ 0.1i? v (where Z ^ 0.1). 

Figure [4] shows the stacked 3D metallicity profile of the 
La emitting gas in the three simulated CDB galaxies. It is 
computed for gas in the temperature range (1—5) x 10 4 K and 
is once weighted by gas density and once by luminosity. The 
metallicity profiles were computed in spherical shells. The 
profile shows that the metallicity falls below 0.1 solar outside 
the disc radius of ~ 8 kpc, and it drops to much smaller 
values outside the inner ~ 30 kpc, where the streams are 
basically made of compressed intergalactic gas. Therefore 
dust opacity is estimated to reduce the La luminosity 
from the streams by a small factor only. However, this 
assumption is likely to fail in the inner galactic disk, 
where the metallicity is above 0.1 solar despite the high 
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Figure 5. Cumulative, mass-weighted temperature distribution 
in the stacked 3 CDB haloes at z = 2.3. The curves refer to 
the streams at r = (0.2 — 1)-R V (dashed blue) and to the whole 
halo including the inner galaxy (solid red). We see that the 
temperature distribution throughout the halo is bimodal, with 
cold and hot phases at 10 4 < T < 3 X 10 4 K and 3 X 10 5 < T < 
2 X 10 6 K, respectively. Outside the inner 0.2_R V , the hot medium 
fills the halo while the cold gas is confined to inflowing narrow 
streams. The cold tail at T < 10 4 K comes from dense clumps 
within the streams. The box marks the temperature range of the 
La -emitting gas. 



redshifts. The effective transmission parameter f a in the 
galaxy vicinity is lower than in the streams. By ignoring 
the dust, we overestimate the La luminosity from the disk. 



5 THE LYMAN-ALPHA EMITTING GAS 

Figure [S] shows the cumulative mass-weighted temperature 
distribution in the CDB galaxies at z = 2.3. The 
distribution is specified alternatively for the halo gas in the 
radius range (0.2 — 1.0)iZ v and for the whole gas, including 
the inner part that involves the disc and its neighbourhood. 
We see that the temperature distribution in the halo is 
bimodal, with a hot virial phase at 10 5 < T < 2 X 10 6 K 
containing ~ 35% of the gas and a cold phase (marked by 
a box) at 10 4 < T < 3 x 10 4 K containing ~ 50% of the 
gas. A very cold tail at T < 10 4 K contains ~ 15% of the 
gas. For the somewhat more massive MN galaxies, which 
we do not show here, the situation is qualitatively similar. 
The temperature distribution is again bimodal, with a hot 
phase at 6 x 10 4 < T < 5 x 10 6 K containing ~ 45% of 
the gas and a cold phase at 2 x 10 4 < T < 4 x 10 4 K 
containing ~ 35% of the gas. An intermediate component 
in the range 4 x 10 4 < T < 6 x 10 4 K contains ~ 20% of the 
gas. In the two kinds of simulations, while the hot phase is 
spread throughout the halo, the cold phase is in the narrow 
streams flowing through the hot medium and the very cold 
tail is concentrated in the dense clumps within the streams 



Figure 6. Cumulative, mass-weighted gas density distribution in 
the stacked 3 CDB haloes at z = 2.3. The curves refer to the 
streams at r = (0.2 — l)R v (dashed blue) and to the whole halo 
including the inner galaxy (solid red). The density in the cold 
streams is in the range n = 0.001 — 0.1 cm -3 , and in the denser 
colder clumps it becomes n = 0.1 — 1cm -3 . The box marks the 
density range of the La -emitting gas, based on the box marked in 
Fig. |51 and assuming that n is a monotonically decreasing function 
of T. 



and in the disc. Based on the temperature dependence of 
the emissivity in Fig. [3] we conclude that most of the La 
luminosity is expected to come from the cold streams, with 
about half the total gas mass participating in efficient La 
emission. 

Figure |S] shows the cumulative mass-weighted gas 
density distribution in the CDB galaxies at z = 2.33 in 
the same radius zones. If the multi-phase medium is in 
pressure equilibrium, the rank order of the cells by density 
and by temperature are exactly opposite of one another, so a 
comparison of Fig.[5]and Fig.[6]reveals that the La -emitting 
cold streams have densities in the range from below 0.01 to 
above 0.1cm -3 (marked by a box). The situation for the 
MN galaxies is comparable. 

Figure [7] shows the two-dimensional gas distribution in 
the temperature-density plane. It complements Figures [S] 
and [6] with information on the different gas phases present 
in the simulations. We see that much of the gas is not in 
pressure equilibrium, nT = const. The top and bottom 
panels display the distribution of mass and La luminosity, 
respectively. We see that the luminosity peaks at T = 
(1 - 2) x 10 4 K over a broad range of densities from n = 
0.001cm" 3 and up. We interpret the T > 10 4 K gas with 
densities above n = 0.1 cm~ 3 outside the central galaxy to 
be in clumps along the streams. It seems to be associated 
with a few percent of the mass that contribute about 25% 
of the luminosity. A similar fraction of the luminosity comes 
form gas in the temperature range T — (2 — 7) x 10 4 K. 

In order to evaluate the possible error introduced in the 
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Figure 7. Distribution of mass (top) and luminosity (bottom) in the n — T plane for the stacked 3 CDB haloes at z = 2.3. The panels 
refer to the streams at r = (0.2 — 1)R V (right) and to the whole halo including the inner galaxy (left). The units are mass Mq (top) 
and erg s _1 (bottom) per unit area in the (logio n, logio T) plane. The luminosity peaks at T = (1 — 2) X 10 4 K over a broad range of 
densities from n = 0.001 cm -3 and up. Most of the emission is from the streams outside the central galaxy. 



MN galaxies by practically ignoring the emissivity from cells 
with densities n > 0.1 cm -3 , we computed the fraction of the 
luminosity coming from such cells in the CDB galaxies. We 
find that this fraction is typically only ~ 15% in the halo at 
r > 0.2R V , and it could be as high as 40% in the central disc. 
By adding 15% to the MN emissivities everywhere we obtain 
a good match between the total luminosities at a given halo 
mass in the CDB and MN simulations (see Fig. I12p . 



6 SIMULATED LYMAN-ALPHA BLOBS 

Figures [8] and [9] show sample images of simulated galaxies, 
two CDB galaxies and two MN galaxies. The haloes are of 



virial mass M v ~ 4 x 10 and 10 Mq respectively, and the 
corresponding redshifts are z — 2.3 and 2.5. 

The top panels present the neutral hydrogen column 
density as computed in ^assuming CIE. The middle panels 
are maps of La restframe surface brightness S, namely a 
fraction f a of the emissivity integrated along the line of 
sight, as emitted per unit area at the galaxy. The bottom 
panels show images of "observed" surface brightness / as an 
observer would see it, per unit area in the galaxy and at the 
telescope. The restframe surface brightness is converted to 
observed surface brightness via 



S 



(20) 



4tt(1 + z) 4 ' 

In order to obtain realistic images, we applied a Gaussian 



© 2009 RAS, MNRAS Q00.[TH2T1 



10 Goerdt, Dekel, Sternberg, Ceverino, Teyssier & Primack 



PSF with a 0.6 arcsec FWHM to mimic atmospheric 
distortions in good seeing conditions, and assumed a pixel 
size of 0.2 arcsec at the telescope. 

One can read from the middle maps the fraction of 
the luminosity that comes from each of the different parts 
of the halo. In particular, we consider separately (a) the 
contribution from the disc galaxy and its near vicinity inside 
a circle of radius 0.2R V and (b) the contribution from the 
halo at 0.2 — 1.0i? v - We see in the figure that the restframe 
surface brightness in the inner galaxy reaches values higher 
than S ~ 10 40 ergs -1 kpc -2 but over a limited area of 
~ 10 2 kpc 2 , thus contributing L ~ a few xlO 42 erg s . 
The typical surface brightnesses in the halo are lower by a 
factor of a few, but the emission regions are spread over an 
area that is larger by an order of magnitude, thus providing 
a comparable contribution to the total luminosity. Most of 
the luminosity comes from the low-density streams, and a 
non-negligible fraction is from clumps associated with the 
streams. 

Recall that the emission from the disc may be partly 
absorbed by dust, which is ignored here, but we do not 
expect substantial dust absorption in the streams in the 
outer halo, where the metallicity is 0.01 to 0.1 solar. We 
can thus consider the halo contribution to be a safe lower 
limit to the overall La luminosity, and take the luminosity 
as estimated from the whole system of disc plus halo as an 
upper limit. 

In the bottom panels we see blobs with surface 
brightness above 10 -18 ergs -1 cm -2 arcsec -2 extending to 
~ 100 kpc. The images are clumpy and irregular. Their 
shapes are non-circular; they tend to be elongated and 
asymmetric, showing finger-like extensions. These sample 
images resemble the ob served images o f LAB s, e.g., those 
shown in Figure 8 of iMatsuda etaH (|2004 ). Figure 1101 
compares "observed" surface brightness maps of two 
simulated CDB galaxies (from the bottom panels of Figures 
[8] and [9]) with t he images of t h e tw o most luminous 
LABs observed by IMatsuda et al.1 (120041 ). The qualitative 
similarity of the irregular, clumpy and elongated LAB 
morphologies and their sizes is encouraging. 



where A, B and C are free parameters, L43 = 
I/Lq/10 43 ergs -1 and M 12 = M V /1O 12 M . To determine the 
free parameters, we first address the mass dependence at a 
fixed redshift where we have performed our most detailed 
analysis in CDB , z ~ 2.4. We then assume that the same 
mass dependence is valid between z ~ 2 and 2 ~ 4, and 
evaluate the redshift dependence. 

Figure [12] shows the luminosity as a function of halo 
mass at z = 3.1. The luminosity is computed within the 
isophotal surface brightness thr eshold of 2.2 x 10 -1 8 erg 
s -1 cm -2 arcsec -2 , also used by IMatsuda et al.1 (|2004T l. We 
adopted f a = 0.85. A small correction upward of ~ 30% has 
been applied to the results as extracted at z ~ 2.4 in order 
to bring them to z — 3.1 in this plot, using eq. (|2ip with 
C = 1.3 as determined below. The average luminosity for the 
three CDB galaxies is shown at M v ~ 4 x 10 11 M o . Each 
of the other symbols is the average of 12 MN galaxies. We 
see that in the mass range most relevant for LABs, M v = 
10 11 ' 5 — 10 12,5 Mq, the luminosity is roughly proportional 
to halo mass. This dependence is driven by the near-linear 
mass dependence of the accretion rate, eq. J5). 

For a more accurate description of the mass dependence 
we fit the function of eq. (|21|l to the results at the fixed 
redshift z ~ 2.4, using least-squares fit in the log, and obtain 
the lines shown in Fig. [12] with B = 0.80 for the whole halo 
and B = 0.76 for the streams only. 

We next use this functional mass dependence to scale 
the simulated results at other redshifts to a fixed mass, 
M v = 1O 12 M . The results are the values as quoted in 
the (M v ,z) bins shown in Fig. 1111 The scaled luminosities 
are shown in Fig. 1131 A crude power-law fit to the symbols 
yields a redshift-dependence power index C = 1.3, with 
the normalisation parameters A — 0.188 and 0.0972 for 
streams+galaxy and for streams only, respectively. We note 
that the obtained redshift dependence of the luminosity is 
somewhat weaker than the (1 + z) 2 " 2 dependence of the 
accretion rate in eq. (0. The accuracy of this scaling with 
redshift is not an important issue for us here since we only 
use it to correct the luminosities from z ~ 2.4 to z — 3.1, a 
correction of less than 30%. 



7 LYMAN-ALPHA LUMINOSITY AS A 
FUNCTION OF MASS AND REDSHIFT 

The variation of average total La luminosity from within 
the virial sphere as a function of halo mass and redshift, as 
derived from the simulations, is shown in Fig. 1111 The two 
panels refer to the luminosity from the whole halo, streams 
and central galaxy, as an upper limit, and to the streams 
at r > 0.2i? v only as a lower limit. The fraction of the 
radiation emitted is set to f a = 0.85, as determined by a 
fit to the observed luminosity function in Ej8] below. The 
luminosities are quoted for small bins about points in the 
(M v ,z) plane. Four points refer to the average luminosities 
over the three CDB galaxies at z — 1.9,2.5,3.2 and 4.0. 
Each of the other points refer to the average over 12 MN 
galaxies. The two-dimensional functional fit to the quoted 
values, L(M V , z), is shown in colours. We use a function of 
the form 

L 43 (M V ) = A {M 12 ) B (1 + zf , (21) 



8 LYMAN-ALPHA LUMINOSITY FUNCTION 

To obtain a predicted LAB luminosity function at z ~ 3, 
we convo lve Z/La(M v ) with the Sheth-Tormen halo mass 
function (|Sheth fc Tormenll 19991 ) at the same redshift. Here 
we assume a cosmology with Qa = 0.72, Qm = 0.28, 
a 8 = 0.83, h = 0.71 and tt B = 0.045. The resultant 
luminosity function is presented in Fig. 1141 both for the 
total luminosity within the virial radius and for the halo 
streams only, in the range 0.2 — 1.0R V , which we consider to 
be upper and lower limits respectively given the uncertainty 
in the dust obscuration from the inner galaxy. The fraction 
of emitted radiation, which could range in principle between 
0.5 and 1.0, is set to f a = 0.85. 

We computed a crude estimate of the observed La 
luminosity function at 2 ~ 3.1 in 3 luminosity bins 
using the most recent 201 LABs observed by Matsuda et 
al. (private communication). This sample contains three 
fields of different environment densities, and can thus be 
crudely considered as a fair sample. The number densities 
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Figure 8. Images of two simulated CDB galaxies at z = 2.3 (left and right). The virial masses are M v ~ 4 x IO^Mq. The box side 
is 140 kpc (physical). The outer circle marks the virial radius and the inner circle is at 0.2R V . Top: Neutral-hydrogen column density. 
Contours are shown for 10 20 and 10 21 cm - 2 . Middle: Restframe surface brightness S, with f a = 0.85, showing contours at 10 39 and 
10 4<) crgs _1 kpc -2 . Bottom: "Observed" surface brightness /, at an angular resolution of ~ 0.1 arcscc, with contours at 10 — 18 and 
10 — 17 ergs _1 cm~ 2 arcsec~ 2 , corresponding to the contours of S. The fraction of luminosity that originates from within these contours 
is 80% and 20%, respectively. 
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Figure 9. Images of two simulated MN galaxies at z = 2.5 (left and right). The virial masses are M v ~ 1O 12 M0. The box side is 
184 kpc (physical). The outer circle marks the virial radius and the inner circle is at 0.2R V . Top: Neutral- hydrogen column density. 
Contours are shown for 10 20 and 10 21 cm -2 . Middle: Restframe surface brightness S, with f a = 0.85, showing contours at 10 39 and 
10 40 ergs — 1 kpc - 2 . Bottom: "Observed" surface brightness /, at an angular resolution of ~ 0.1 arcsec, with contours at 10 -18 and 
10 — 17 ergs -1 cm -2 arcsec -2 , corresponding to the contours of S. The fraction of luminosity that originates from within these contours 
is 93% and 18%, respectively. 
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Figure 10. A qualitative crude comparison of im ages of simulated (top) and observed (bottom) LABs. The observed LABs are the two 
most luminous LABs from Matsuda ct al. (2004) with L ~ 10 44 crgs -1 . The simulated images are of CDB galaxies scaled from z = 2.3 
to z = 3.1 and to the luminosity of the corresponding observed LABs using I oc Roc L 1 / 3 (the original, unsealed luminosity was 5 X 10 43 
erg s _1 ). All four maps have an angular resolution of ~ 0.2 arsec. 



were determined by counting the number of blobs in 
the luminosity bins and dividing by the comoving survey 
volume of V s = 1.5 X 10 6 Mpc 3 (or physical volume 
3.6 x 10 5 Mpc 3 ). Also shown is the luminosity function 
derived earlier from a partial sample o f 35 L ABs in a dense 
cluster environment (|Matsuda et al.1 l2004l i 3 l. The latter 
overestimates the universal luminosity function by a factor 
of a few as expected from the high environment density. A 
Poissonian error has been attached to every bin. The number 



3 The comoving survey volume here is V s = 1.3 X 10 5 Mpc 3 (or 
physical volume 3.2 X 10 4 Mpc 3 ). They used a cosmology with 
Om = 0.3, £7a = 0.7 and h = 0.7, which is very close to the 
cosmologies used in the simulations and in the analysis. They 
choose their sample as explained in their section 3. The observed 
luminosity function was computed for a combined limit of surface 
brightness (2.2 X 10~ 18 erg s" 1 arcsec -2 ) and size (16 arcsec -2 ), 
while in the simulations we use only a surface brightness limit. 
Had we imposed a similar size limit in the simulations, one 
can deduce from Fig. 1161 that our predicted luminosity function 
would have been suppressed at luminosities below a few times 
10 42 ergs -1 , making the comparison with the counts at the faint 
end less certain. 



density has been derived by dividing the number of objects 
in each bin by the total survey volume. 

We learn from Fig. [14] that the predicted luminosity 
function has a similar slope to the observed one. For 
f a ~ 0.85, the predicted upper and lower limits border 
the observed function from above and below throughout the 
whole range of LAB luminosities. We find that f a changes 
roughly linearly with the adopted cosmological erg through 
its effect on n(M). Therefore, the ~ 10% uncertainty in erg 
translates into a similar uncertainty in f a . 

Also shown is the luminosity function as estimated from 
our crude toy model for gravitational heating, ij2] again 
convolved with the Sheth-Tormen halo mass function, and 
using fafc = 0.34. It is remarkable that this very crude toy 
model recovers the simulation results in the relevant mass 
range to within a factor of a few. 



9 SURFACE-BRIGHTNESS PROFILE, AREA 
AND LINEWIDTH 

A detailed comparison of theory to data will be presented 
in a future paper where we apply a more accurate radiative 
transfer calculation to the simulated galaxies. However, it is 
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Figure 11. Total La luminosity as a function of halo mass and redshift, separately for streams+galaxy (left) and for streams only 
(right). The average values as drawn from the simulations in bins of (M v , z) are marked by symbols and numbers in units of 10 43 ergs -1 . 
The CDB results are marked by filled symbols and the MN results by open symbols. Shown in colour is the fitting function eq. ( 12 U with 
contours at 10 42 , 10 43 and 10 44 erg s" 1 . 



worthwhile to report here on a preliminary comparison with 
data, based on the simplified analysis of the current paper, 
which is encouraging. 

9.1 Surface-brightness profile 

Figure [151 shows the surface-brightness profile of the stacked 
images from the three CDB simulations, assuming f a — 0.85 
throughout the whole halo. It is compared to stack ed profiles 
of 35 observed LABs from iMatsuda et all (|2004 l. For the 
stacking of the observed LABs, which range over more than 
an order of magnitude in luminosity, we scale the LAB radius 
R and surface brightness 7 to a fiducial luminosity L = 10 43 
erg s . We write L oc I R 2 and assume L oc M v and R oc 
R v oc M, 1 ' 3 to obtain the scaling R oc I oc L 1//3 . Thus, for a 
LAB with luminosity L, we multiply both the radius and the 
surface-brightness by the factor (Lo/L) 1 ^ 3 . Non-zero values 
for I are considered for each LAB only above the isophotal 
surface-brightness threshold of 2.2 x 10 -18 erg s _1 cm -2 
arsec -2 and only the inner adjacent non-zero values (i.e. all 
values outside the innermost zero value are also treated as 
zero values) were used in the averaging of I at each radius 
r. The stacked simulated profile has been scaled accordingly 
to match the same fiducial luminosity Lq. 

We see that a power law I(r) oc r -7 is a good fit to the 
simulated profile, with 7 ~ 1.2, from the disk scale r < 10 
kpc out to ~ 40 kpc, which is about half the virial radius. 
It seems to steepen to 7 ~ 2 a t larger radiiQ The observed 
images of lMatsuda et all (|2004h are subject to a point spread 
function of ~ 1 arcsec, which is responsible for the flattening 

4 Note that the slope 7 ~ 1.2 is consistent with the slope of the 
cumulative 3D luminosity profile shown in Fig.[T] where L(< r) oc 
r - 8 in the range r = (0.1 - 0.4) R v . 



of the observed profile at r ^ 10 kpc. Given the high 
background outside the surface-brightness threshold, the 
meaningful part of the observed stacked profile is limited 
to the range 10 — 30 kpc. In this range, the power law with 
a slope 7 = 1.2 provides a good fit to the observed profile as 
well. We also see that the crude toy model of ^provides a 
profile that is consistent with the simulated profile to within 
a factor of two, as seen in Figure [T] 



9.2 Isophotal area 

Two of the observable global quantities for each LAB 
are the area encompassed by an isophotal contour of a 
given threshold surface brightness, and the corresponding 
lumin osity. In order to compare to the data of lMatsuda et all 
(|2004r ), we applied to our simulated galaxies, scaled to 
z = 3.1, a threshold of 2.2 x 10 -18 erg s _1 cm -2 arcsec -2 . 

Figure [TS] displays the areas and luminosities of our 
simulated LABs in comparison with the 35 observed LABs. 
Each of the three simulated CDB galaxies is "observed" 
from three orthogonal directions, and each of the randomly 
selected MN galaxies is "observed" from one random 
direction. We see that the area at a given luminosity in 
MN is systematically smaller than in CDB , but only by a 
factor of two. The simulations agree with the observations 
to within a factor of two, and they reproduce the general 
correlation between area and luminosity. The overall trend 
follows the scaling relation 

A oc L 0,75 . (22) 
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Figure 12. La luminosity as a function of halo mass at z = 
3.1 with / a = 0.85, for the total luminosity within the virial 
radius (upper, red symbols and curve), and for the halo outside 
0.2R V (lower, blue symbols and curve). Each symbol represents 
the average luminosity over simulated haloes of a given mass, 12 
haloes from the MN simulation (open symbols) and 3 from the 
CDB simulations (filled symbols). The lines are least-squares fits 
using eq. I|21j l with slopes B = 0.80 and 0.76, respectively. 



9.3 Linewidth 

The major contribution to the LAB line profile is the 
kinematic effect owing to the distribution of line-of-sight 
velocities in the cold gas. It is dominated by the coherent 
instreaming velocities from the background and from the 
foreground, corresponding to a FWHM of order a few times 
the halo virial velocity. The effect of resonant scattering of 
La should be convolved with the kinematic effect. 

In Fig. [17] we show the kinematic components of three 
characteristic line profiles from the simulated CDB galaxies. 
The effect of resonant scattering is not included yet. For 
any desired direction of view, the kinematic line profile is 
computed as the luminosity weighted distribution of the 
line-of-sight velocities across the whole isophotal area. In 
some cases, the line profile is dominated by one peak, 
corresponding to a stream that stretches roughly along the 
line of sight. In other cases, the line profile is bimodal, 
with a red peak and a blue peak corresponding to two 
opposite streams that lie not far from the line of sight. A 
third type of line profile is flatter, showing three peaks, the 
central of which corresponding to a stream that is roughly 
perpendicular to the line of sight. We crudely estimate the 
kinematic FWHM as twice the standard deviation of An. 

Resonant scattering of the La photons in the stream 
neutral hydrogen gas broadens the La emission line profile 
as the trapped p hotons diffuse into the damp ing wings before 
finally escaping (|Adamsll 19721 lNeufeldlll990| ). For scattering 
and escape from a static plane parallel medium, t he total 
line width, in units of the Doppler parameter, is l|Adamsl 
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Figure 13. ha luminosity as a function redshift for M v = 
10 12 Afo with f a = 0.85, for the total luminosity within the 
virial radius (upper, red symbols and curve), and for the halo 
outside 0.2i? v (lower, blue symbols and curve). Each symbol 
represents the average luminosity over simulated haloes of a given 
mass, 12 haloes from the MN simulation (filled symbols) and 3 
from the CDB simulations (open symbols). The luminosities were 
scaled from the simulated values at the different halo masses to 
M v = 10 12 M Q using the mass dependence obtained in cq. J2H . 
The curves are approximate eye-ball power-law fits with a slope 
C= 1.3. 



2[ar LQ 7(3/^)] 1/3 



(23) 



For La , the Doppler parameter is b — 12.89 T^ 2 km s~ x 
and the radiation parameter a and the optical depth tl 
are given in Sj4] following eq. (|19[) . We obtain 



436A4 /3 T 4 1/6 kms 



(24) 



ll972l ; lHarringtonlll973l ) 



Thus, the line profile resulting from resonant scattering of a 
midplane source is expected to consist of a red and a blue 
peak separated by <T scat . For a very crude estimate of total 
line width An, we add in quadrature <T SC at to the FWHM of 
the kinematic line profile|j 

Figure [18] displays isophotal area versus linewidth 
for the simula t ed L ABs and for observed La emitters 
iMatsuda et~aH (|2006l ). including LABs and less extended 
emitters. The simulations and observations show a similar 
general trend, which very crudely follows A oc (AV) 2 ' 25 , 
but the data show large scatter. A comparison to the 
tighter correlation seen in Fig. [16] between the area and the 
luminosity indicates that the scatter in Fig. [18] is mostly 
due to scatter in Av. This scatter is partly due to the 
variations in the relative orientations of the streams and 
the line of sight. The spread in the Av as estimated in the 
simulations is underestimated due to the fact that we added 
in quadrature a constant value of <T sca t = 436kms _1 . This 

5 We note that the effect of resonant scattering as determined 
for a static medium may be an overestimate because of turbulent 
motions. 
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Figure 14. ha luminosity function, showing the comoving 
number density of LABs brighter than L^ a at z = 3.1. The 
observed symbols (black, filled circles) are based on a sample 
of 201 LABs by Matsuda (private communication) in a fair 
sample. The earlier, higher estimates (magenta, open circles) 
are based on 35 LABs in a dense cluster environment. The 
lower and upper limits for the luminosity function as derived 
from the simulations are shown, for streams only (dashed, blue 
curve) and for strcams+galaxy (solid, red curve), respectively. 
The transmission factor was set to f a = 0.85 for a best 
match between theory and observation. The gravitational-heating 
toy-model prediction based on cq. @ is shown (green, dashed 
line). 

is also responsible for the artificially sharp lower bound at 
Av = 436 km s" 1 , which is apparent for small values of A. 

Finally, Fig. [19] displays linewidth versus luminosity 
for the simulations and a subset of the observed LABs 
by iMatsuda et all |2006l ). The simulations reveal a scaling 
relation that can be approximated by Av oc L 018 . The 
simulations and observations agree to within a factor better 
than two, but the larger scatter in the observed data 
confirms our suspicion that we have underestimated the 
scatter in the analysis of Av from the simulations. 

We conclude that our simplified analysis of La 
emission from simulated galaxies qualitatively reproduces 
the observed correlations between the global quantities of Lev 
luminosity, area and linewidth. A more detailed comparison 
of theory and observed LABs, especially involving line 
profiles and linewidths, should await a more accurate 
analysis of the simulations using radiative transfer and 
including dust absorption, as well as photometric and 
spectroscopic measurements of a larger sample of observed 
LABs. 



10 GRAVITATIONAL HEATING 

We now return to the energy source for the La emission 
in our simulations. In addition to the gravitational 
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Figure 15. ha surface-brightness profile. The stacked profile 
from the three CDB galaxies (red, solid line) scaled from z = 2.3 
to z = 3.1 is compared to t he stacked profile of 35 observed LABs 
from lMatsuda et alj i20pj) at z = 3.1 (blue, dotted curve). The 
individual profiles are scaled to a LAB of a fiducial luminosity of 
10 43 erg s — 1 (see text). The toy- model prediction from section[2]is 
shown (green, dashed curve). A power law of slope —1.2 is shown 
as a reference (pink, dotted line). There is a good agreement 
between simulation and data in the radius range 10 — 30 kpc, 
where the observed profile is meaningful. 




Figure 16. Isophotal area versus luminosity of LABs at z = 3.1, 
above a surface-brightness threshold of 2.2 X 10 — 18 erg s — 1 cm~ 2 
arcsec - 2 . The simulated LABs from CDB (red filled circles) and 
MN ( green open circles) are compared with the 35 observed LABs 
from IMatsuda et al. | | |2004) (blue filled squares). The pink solid 
line refers to a power law of slope 0.75, which resembles the trend 
seen both in the simulations and the observations. 
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Figure 17. The kinematic components of three characteristic La 
line profiles from CDB simulated galaxies, tentatively neglecting 
the effect of resonant scattering. These examples have between 
one to three dominant peaks, reelecting the relative orientations 
of the big cold streams and the line of sight. 
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Figure 18. Isophotal area versus lincwidth (FWHM) for LABs 
at z = 3.1. The simulated LABs from CDB (red filled circles) 
and MN (green open c ircles) are compared with La emitters from 
iMatsuda et alj ll2006h (blue filled squares). The pink solid line 
refers to a power law of slope 2.25, which crudely resembles the 
general trend seen in the simulations and in the observations. 
The Aw deduced from the simulations is an underestimate and is 
subject to an artificially sharp lower bound at 436 km s . 
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Figure 19. La luminosity versus linewidth (FWHM) for LABs 
at z = 3.1. The simulated LABs from CDB (red filled circles) 
and MN (green open circle s) are compare d with LABs from 
IMatsuda et al.l (120041 1 and IMatsuda et al.l j2006l l (blue filled 
squares). The pink solid line refers to a power law of slope 0.18, 
which crudely resembles the general trend seen in the simulations 
and in the observations, which are consistent within a factor of 
two. 



heating, which was crudely approximated in [J2] there is 
photoionisation by the UV background. Photoionisation 
effects by central sources such as stars and AGN were not 
simulated, and they are not expected to be substantial 
because of the shielding implied by the radial orientation 
of the cold streams. Before we address gravitational 
heating again, we should verify the contribution of the UV 
background. 

A simple estimate is as follows. For an isotropic 
background with mean photon intensity 47tJ* and mean 
photon energy e, the rate at which energy is absorbed 
across a spherical surface of radius R is 7v,re4TrR 2 . If we 
adopt for the metagalactic field at z = 3 a flux ttJ* ~ 
5.5 x 10 4 photons s _1 cm -2 and e ~ 20eV, we obtain a total 
UV heating rate into the virial radius R v ~ 100 kpc of 
2 x 10 42 ergs~ 1 . Even if all of the UV energy goes into 
heating the cold streams, the UV heating is small compared 
to the total La luminosity of ~ 10 43 ergs _1 and to the 
gravitational heating rate as estimated in However, 
one should keep in mind that the MN simulation may 
overpredict the energy coming from the UV background 
when assuming that the gas is optically thin and the UV 
background is uniform. 

To evaluate the relative contributions of gravitational 
heating and UV background in our actual simulation, we 
read from each grid cell of the simulation snapshot the 
instantaneous input rate of energy per unit volume by 
the UV flux, euv, and subtract it from the La emissivity 
computed in eq. JSJ, to obtain the relative contribution of 
gravitational heating in that cell, 



© 2009 RAS, MNRAS 000,rnf2T1 



18 Goerdt, Dekel, Sternberg, Ceverino, Teyssier & Primack 



0.8 



0.6 



0.4 



0.2 



streams + galaxy 




streams 



0.2 0.4 0.6 0.; 

fgrav 



Figure 20. The role of gravitational heating versus 
photoionisation. Shown is the cumulative, luminosity-weighted 
distribution of the fraction of the energy that is provided by 
gravitational heating, /grav, in the three CDB galaxies. The 
curves refer to the streams at r = (0.2 — l)R v (dashed blue) and 
to the whole halo including the inner galaxy (solid red). The 
area under the curve implies that for the total CDB halo within 
Rv 86.8% of the energy input is gravitational, and only 13.2% is 
due to the UV background. The numbers are 84.2% and 15.8% 
respectively for the CDB halo between 0.2 and 1.0 -R v - 



euv 



(25) 



Figure [20] shows the cumulative luminosity-weighted 
distribution of / gra v for the three CDB galaxies at z = 2.3. 
In this plot /grav = refers to cells where all the La 
luminosity is driven by the UV background and /grav = 1 
corresponds to gravitational heating only. The area under 
the curve is the fraction of the total energy provided by 
gravity. In the CDB haloes this fraction is 86.8%, with only 
13.2% due to the UV background for the total halo within 
R v . It is 84.2% and 15.8% respectively for the halo between 
0.2 and 1.0 R v . In more detail, if we focus, for example, 
on the cells of gas with the highest values of / gra v that 
contribute 0.8 of the total luminosity, we read from the 
figure that they all have /grav > 0.8. This means that in the 
gas that is responsible for 80% of the La luminosity, more 
than 80% of the energy is gravitational and less than 20% is 
due to the UV background. The fraction of the total energy 
provided by gravity in the MN haloes is 91.4%, with only 
8.6% due to the UV background for the total halo within 
R v . It is 79.5% and 20.5% respectively for the halo between 
0.2 and 1.0 R v . We conclude that most of the La emission 
from our simulated galaxies is indeed driven by gravitational 
heating. 

The total luminosity as predicted by the simplified 
gravitational-heating toy model can be compared to the 
total La luminosity that is actually produced in the 
simulations. We focus on the halo streams in the radius 



range r = (0.2 - 1.0) R v in haloes of M v = 1O 12 M at 
z = 3. The toy-model estimate, based on eq. (J7J) and Fig. [T] 
is -Ehoat ~ 10 3 ergs -1 . A straightforward integral of the 
emissivity in the simulated galaxies over the same volume 
in the halo yields a comparable value (which can be read 
as L/ f a from Fig. I12p . This implies that the toy-model 
parameter / c should be of order unity, indicating that most 
of the gravitational energy is deposited in the cold streams. 

We now return to Fig. [T] where we showed 
the gravitational heating toy-model prediction for the 
cumulative 3D luminosity profile. Shown in comparison is 
the corresponding La luminosity profile as derived from 
the actual simulated CDB galaxies, averaged over the three 
galaxies. The similarity of the two profiles is remarkable. 
This implies, again, that the La emission is indeed powered 
by gravitational heating, where the potential energy of the 
instreaming cold gas is radiated as La . 

The cumulative luminosity profile can be fitted by a 
power law, 



L(< r) oc r° 



(26) 



corresponding to a 3D luminosity density profile £(r) oc 
r a ~ 3 . In the range r = (0.1 — 0.4)i? v , we find for the 
simulated luminosity a ~ 0.8, namely £ oc r~ 2 ' 2 . This is 
consistent with the average surface-brightness profile I(r) ~ 
r -1 ' 2 measured in H9.ll The toy model predicts a similar 
power index over the same radius range. We also note that 
the density profile of cold gas follows a similar power law. 



11 DISCUSSION AND CONCLUSION 

Hydrodynamical cosmological simulations 

robustly demonstrate that massive galaxies at high redshifts 
were fed by cold gas streams, inflowing into dark-matter 
haloes at high rates along the cosmic web (iKeres et al.ll2005l ; 
iDekel fc Birnboiml 120061 ; iDekel et al.1 120091 ). In this paper 
we have shown that these streams should be observable as 
luminous La sources, with elongated irregular structures 
stretching for distances of over 100 kpc. The release of 
gravitational potential energy by the instreaming gas as it 
falls into the halo potential well is the origin of the La 
luminosity, which ranges from 10 42 to 10 44 ergs -1 for haloes 
in the mass range 10 11 - 10 13 M Q at z ~ 3. The predicted La 
emission morphologies and luminosities make such streams 
likely candidates for the sources of observed high redshift La 
blobs. Most of the gas in the cold streams is at temperatures 
in the range (1 — 5) x 10 4 K, for which the La emissivity is 
maximised. The hydrogen gas densities in the streams are in 
the range 0.01 — 0.1 cm 3 , and are higher in clumps that flow 
with the streams, leading to the high surface brightnesses. 

Using state-of-the-art cosmological AMR simulations 
with 70 pc resolution and cooling to below 10 4 K, and 
applying a straightforward analysis of La emissivity due 
to electron impact excitation, we produced maps of La 
emission from simulated massive galaxies at z ~ 3. We 
computed the average La luminosity per given halo mass 
and the predicted luminosity function of these extended La 
sources, with an uncertainty of ±0.4 dex. The properties 
of the individual images resemble those of the observed 
LABs in terms of morphology and kinematics. The predicted 
luminosity function is close to the observed LAB luminosity 
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function. The simulated LABs qualitatively reproduce the 
observed correlations between the global LAB properties of 
La luminosity, isophotal area and linewidth. 

The LAB properties can be understood using a very 
simple toy model that accounts for gravitational heating of 
the inflowing gas. In comparison, the more elaborate toy 
model of DL09 predicts a steeper power-law relation between 
luminosity and halo mass (their Fig. 3, with / g = 0.2). In 
the mass range 10 12 — 1O 13 M0, DL09 predict luminosities in 
the range 1.5 x 10 42 — 2.5 x 10 44 ergs -1 . A comparison to our 
Fig. [12] indicates that they underestimate the luminosity at 
M v ~ 10 12 M Q by a factor of a few. The DL09 toy model 
matches the observed luminosity function only after the 
predictions are boosted, trying to account for an overdense 
region of the universe. The association of the predicted 
LABs with massive haloes places the m preferentially in 
overdense environm e nts, a s observed (|Steidel et al.l |2000| ; 
Mat suda et al. 1 12004 I 20061 ). but this is already accounted 
for by the number density of haloes in a fair sample. 

The association of observe d LABs with star-forming 
LBGs jHavashino et all |2004). sub- m illimetre galaxies 
JChapman et alj|200ll; iGeach et al.ll2005l) or active g alactic 



nuclei ( Geach et al. 12009^ iBasu-Zvch fc ScharjlioM ) is not 
surprising since star formation and AGN activity are 
triggered by the same process of streaming of cold gas 
into massive haloes. La emission can be driven by any of 
these sources of energy. For example, in photoionised gas 
in star-forming regions the La luminosity is Ll q = 1-5 x 
10 42 SFRergs -1 where SFR is the star-formation rate (in 
Mq yr -1 ) for co ntinuous star formation for a Kroup a initia l 
mass function |Sternberg. Hoffmann fc Pauldrachl 12003? ). 
Thus the typical LAB luminosity of 10 43 erg s _1 would 
require SFR ~ 10 Mq yr -1 . However, most of this radiation 
is likely absorbed by dust, and confined to the central galaxy. 

A limitation of our AMR simulations arises from the 
artificial pressure floor imposed in order to properly resolve 
the Jeans mass. This may have an effect on the temperature 
and density of the cold gas in the streams, with potential 
implications on the computed La emission. Still, the AMR 
code is the best available tool for recovering the stream 
properties. With 70 pc resolution and proper cooling below 
10 4 K, the CDB simulations provide the most reliable 
description of the cold streams so far. The rather small 
correction that we had to apply to the luminosities extracted 
from the MN simulation indicates that the MN galaxies can 
be used to recover the mass and redshift dependence of the 
global stream properties and the scaling relations between 
them. 

Another source of uncertainty has to do with the 
simplified way the ionisation by the UV background is 
handled in the simulations, and with the post-processing 
calculation of the ionisation state of the cold gas. In the 
CDB simulations, the centres of the streams, where the 
gas density is higher than 0.1cm -3 , were assumed to be 
self-shielded against the UV background, in agreement 
with our analytic estimates in ij4] The ionisation state of 
the gas, which is a key ingredient in evaluating the La 
emissivity, was computed in post-processing assuming CIE. 
Cooling rates were computed for the given gas density, 
temperature, metallicity, and UV background based on 
Cloudy (CDB ) or assuming ionisation equilibrium for H 
and He, including both collisional- and photo-ionisation 



(MN). Photoionisation is neglected since the streams 
are sufficiently thick to be self-shielded. An alternative 
calculation using Cloudy yielded lower fractions of neutral 
hydrogen at n ^ 0.01cm -3 , and total luminosities that are 
typically three times smaller than obtained using CIE. This 
calculation assumed that each cell, with a given hydrogen 
density n, is in the middle of a uniform slab of thickness 1 
kpc. Based on our estimates of self-shielding, eq. |T7J), we 
adopt the CIE results as the more reliable approximation. 

In our computations so far we did not consider the 
radiative transfer of La photons through the streams, or 
through the intervening intergalactic medium. We assume 
that most of the radiation emitted at z ~ 3 will reach 
the observer at z — without undergoing significant 
attenuation. Resonant line scattering in the optically thick 
streams will tend to spread out the La emission region, while 
the repeated Doppler shifts broaden the line profiles, and 
the increased path length of the random walk amplifies the 
absorption by dust. Such effects are expected to modify the 
images and line profiles, but to have only a small effect on 
the total La luminosity. An analysis of La emission from our 
simulated galaxies including the effects of radiative transfer 
and dust absorption and a more accurate treatment of 
photoionisation from stars will be presented in a forthcoming 
paper (Kasen et al, in preparation). 

Other hydrogen emission lines, such as Lf3 or Ha, are 
expected to be tw o orders of magnitude less luminous than 
La in our model (|Baldwinlll977l ; lMilledll974h . The column 
densities of Ahi ~ 10 20 cm -2 in the cold streams should also 
be de t ectable as La absorption in quasar s pectra (|Prochaskal 
1 19991 ; I Wolfe. Gawiser fc Prochaskal I2005D . Considering all 
the galaxies that are intersected by a line of sight to a 
background quasar with an impact parameter < R v , we 
crudely estimate from the simulations that absorption by 
N20 > 1 hydrogen should be detected in ~ 20% of the 
galaxies. Alternatively, La photons that are emitted form 
the central galaxy should be absorbed in the radial streams 
feeding the same galaxy at A20 > 10, but only in ~ 5% 
of the galaxies. The streams coul d also be detectable as 
low-ionisation metal absorbers (e.g. lGnat fc Sternberd2007T ) 
as long as the metallicity in the streams is greater than 
~ 0.01 solar (paper in preparation). 

Our results support the idea that the observed LABs 
are direct detections of the cold steams that drive the 
evolution of massive galaxies at high redshifts. Even though 
the observed LABs are sometimes associated with central 
sources that are energetic enough to power the observed 
La emission, such as starbursts and AGNs, these central 
sources are very different from each other in the different 
galaxies, and in many LABs they a re absent altogether 
(|Yang et al] 120081 ; IGeach et all 120071 ). The gravitational 
heating associated with the inflowing cold streams is a 
natural mechanism for driving the extended La cooling 
radiation observed as LABs, and this extended La emission 
is inevitable in most high-redshift galaxies. 
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